Preprint typeset in JHEP style - HYPER VERSION 



CERN-PH-TH /2005-222 
hep-ph/0601075 



Hard pion and prompt photon at RHIC, 
from single to double inclusive production 



Francois Arleo 

CERN, PH-TH Dept. 

1211 Geneva 23, Switzerland 



E-mail: arleo@cern.ch 



Abstract: Single pion and prompt photon large transverse momentum spectra 
in p-p and Au-Au collisions are computed in perturbative QCD at RHIC energy, 
^s NN = 200 GeV. Next-to-leading order calculations are discussed and compared 
with p-p scattering data. Subsequently, quenching factors are computed to leading 
order for both pions and photons within the same energy loss model. The good 
agreement with PHENIX preliminary data allows for a lower estimate of the energy 
density reached in central Au-Au collisions, e RHIC > 10 GeV/fm 3 . Double inclusive 
7-7T production in p-p and Au-Au collisions is then addressed. Next-to-leading 
order corrections prove rather small in p-p scattering. In Au-Au collisions, the 
quenching of momentum-correlation spectra is seen to be sensitive to parton energy 
loss processes, which would help to understand how the fragmentation dynamics is 
modified in nuclear collisions at RHIC. 
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1. Introduction 

Photon production is a promising observable in high energy heavy-ion collisions. 
On the one hand, one expects that hot media such as quark-gluon plasma radiate 
thermal photons [1-4] with transverse momenta of the order of its temperature, 
p ± ^ = 0{T). The huge decay background from neutral pions, especially below 
4 GeV, makes the experimental extraction of such a signal quite a difficult task. In 
that respect, the precise PHENIX measurements recently shown at Quark Matter 
down to 1 GeV are remarkable [5,6]. What is more, the possible excess in photon 
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production reported in the 1-3 GeV range above the expected QCD rate is exciting; 
yet any definite conclusion would be highly premature. On the other hand and 
side of the spectrum, the hard prompt photon signal with p ± ^ ^> T can be used 
to gauge the strength of nuclear effects in hard processes (see e.g. Refs. [7-13] on 
prompt photon phenomenology). In particular, comparing blind hard probes - Drell- 
Yan, heavy bosons, prompt photons - with coloured hard probes - hadrons, jets - 
would allow for a clear experimental distinction between initial-state effects (such as 
small-x saturation physics [14]) and final-state interactions (e.g. parton energy loss 
processes [15,16]). However, it may be misleading to assume that prompt photon 
production, because of its colour neutrality, should not depend on the dense-medium 
properties. Indeed, to leading order in the perturbative expansion, prompt photons 
can be produced directly in the hard subprocess ( "Drell-Yan-like" ) but also from 
the collinear fragmentation of a hard quark or gluon ("jet-like") [11]. Of course, 
only the sum of these two components is meaningful and as scale-independent as 
possible: the leading-order (LO) fragmentation process may be seen as a next-to- 
leading order (NLO) direct contribution and vice versa, depending on the resolution 
scale. Nevertheless, it is worth stressing that the prompt photon signal could in 
principle be modified by the dense-medium produced in the nucleus-nucleus reaction. 
In particular, the mechanism that spectacularly quenches the pion p^ spectra in 
Au-Au collisions [17-19] should also affect prompt photon production, although not 
necessarily with a similar strength. It is thus important to treat on an equal footing 
large-p ± pion and prompt photon production, from p-p to A-A collisions. 

Single hadron spectra are certainly useful to reveal the formation of a dense 
medium in heavy-ion collisions (as we shall see, the present data actually allow for 
a lower estimate of the energy density reached in central Au-Au collisions), however 
such measurements hardly inform us on how exactly the medium affects fragmenta- 
tion functions and, more generally, on the fragmentation process itself. Indeed, the 
initial parton momentum k ± , - hence the fragmentation variable z = p ±h /k ± enter- 
ing fragmentation functions - is not fixed. As a consequence, there is a clear need to 
go beyond the single particle production picture. This triggered in particular several 
measurements [20-24] and calculations [25,26] of 2-particle azimuthal correlations. 
Moreover, following Ref. [27], performing photon-tagged momentum correlations in 
the double inclusive production of 7-7?° pairs at the LHC has been proposed as a 
powerful tool to extract (or at least to constrain) fragmentation functions: to leading 
order, the hard photon (hopefully produced directly) gives access to the leading par- 
ton transverse momentum, which eventually fragments into the pion [28]. The wealth 
of RHIC large-p ± data provides hints that the energy loss process is probably at work 
in Au-Au central collisions [29-31] (and to a lesser extent in Cu-Cu [32]), such a 
picture being supported by many phenomenological analyses [33-37]. Consequently, 
addressing such 7-7?° momentum correlations also at RHIC energy is particularly 
relevant. 
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In this paper, we first explore parton energy loss effects on both single hard pion 
and single prompt photon production within the same model, based on the Baier- 
Dokshitzer-Mueller-Peigne-Schiff (BDMPS) framework [38,39]. After discussing 
single spectra in p-p collisions computed in QCD at NLO, we predict the expected 
pion and photon quenching factors and compare them with the available data. The 
photon total yield over the background is investigated as well. A lower estimate 
for the RHIC energy density is then determined. The second part of this work is 
devoted to the study of photon-tagged momentum correlations. Various absolute 
correlation spectra are constructed in p-p and their expected quenching in central 
Au-Au collisions is discussed. Counting rates are also given. Let us now start with 
the details of the perturbative calculations. 



2. Preamble: predictions 

We discuss in this section the accuracy as well as the ingredients used in the present 
perturbative QCD predictions of single and double inclusive pion and photon pro- 
duction in p-p and Au-Au collisions. 

2.1 Proton— proton collisions 
2.1.1 Single inclusive production 

Single inclusive pion and photon hadropro duct ion cross sections in p-p collisions are 
computed at NLO accuracy in QCD, using the work of Ref. [11]. For single-pion 
production, NLO cross sections read [11] 

da* x 



dp±dy . ^ 



X 



J dx 1 dx 2 F i/p (x 1: M)F j/p (x 2: M)^D v/k {z, M F ) 



2n J dp^dy \ 2n 



K ij:k (v,M,M F ) 



(2.1) 



where Fij/ P are the proton parton distribution functions (PDF), D n /k the pion frag- 
mentation functions, and a^ (respectively the leading-order (respectively next- 
to- leading order) partonic cross section in the MS scheme 1 . We denote by fi, M and 
Mp the renormalization, the factorization and the fragmentation scale. 

Unlike pions, single prompt photons can be produced "directly" in the hard 
subprocess [11], with a contribution 

J 0-7 dir r 

= / davits F i/p (x u M) F j/p (x 2 ,M) 



dp.dy 

t ^ ± y i, 3 =q,g ■ 



2tt V d P± d 2/ 27r ij 



+ ^K?n^M,M F )\, (2.2) 



1 We omit the explicit dependence of afj and on the kinematic variables x\, x 2 , ^/s, p ± , 

and y for clarity. 
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in addition to the collinear fragmentation process [11] 
da 7fmg 



i,3,k=q,g 



/d z 
dxidx 2 — Fi/p(x!,M) Fj/ p (x 2 ,M) D l/k (z 1 M F ) 

{Si + ^(-.-4 (2 ' 3) 

at the same order in the perturbative expansion. Let us once more emphasize that 
the distinction between direct and fragmentation photons is arbitrary, only the sum 
of these two contributions being meaningful and with a lesser fragmentation scale 
dependence [11]. 

2.1.2 Double inclusive production 

At leading order in QCD, the basic two-particle 7-71"° correlation cross section, from 
which various observables are constructed, can be written [40,41]: 

d(T PP ^ 1 Ar/c(*7r, M p ) Dj/djZj, M F ) 



dp^dyAzAp^dy^dz^ Ms 2 a " c > z % 

(2.4) 



F a/p {x!,M) F b/p (x 2 ,M) 



Xi x 2 

where Ai is the LO hard ab — > cd scattering amplitude and the distinction between 
the direct and the fragmentation component is made implicit in the D J u(z 1 , M F ) 
fragmentation functions. 

2.2 Nucleus— nucleus collisions 

The status of NLO QCD calculations in nucleus-nucleus collisions is not as yet 
well established. Therefore, we shall consider LO calculations and only show the 
normalized ratio of the Au-Au over the p-p production cross section: 

1 a do** I dtp* 

where <7 geo is the geometric cross section obtained via the Glauber multiple scat- 

AuAu L 

tering theory, <t nn the nucleon-nucleon cross section, and (iV coll )| the number of 
binary collisions in a given centrality class C. Numerically, we have af™ = 6900 mb, 
<r NN = 42 mb and (N coll )\ c = 779 in central (C < 20%) Au-Au collisions at RHIC 
energy [4, 19]. 

Ignoring any nuclear effect, the nucleus-nucleus collision cross section is deduced 
directly from the p-p LO cross section, replacing the proton PDF in Eqs. ( |2.1| ) to 
by 

F i/A (x,M) = ZF i/p (x,M) + (A-Z)F i/n (x,M), (2.6) 
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where Z and A are respectively the number of protons and the atomic mass number 



of each nucleus. The neutron PDF Fi/ n in Eq. ( |2.6|) is obtained from the proton Fi/ P 
by the usual isospin conjugation assumptions: u p = d n , d p = u n , u p = d n , d p = u n , 
and s p = s n . 

Such possible isospin effects, when comparing different nuclear targets - and 
in particular almost isoscalar nuclei (A ~ 2 Z) such as Au with a proton (A = 
Z) - may be significant when hadron or photon production occurs at large Bjorken 
x = 0(2p ± /yf~s) < 1, at which the partonic process involves essentially the scat- 
tering of valence quarks. Conversely, it should remain completely negligible, say, 
around x ~ 0.01, below which the nucleon PDF is dominated by the gluons. Note 
that such an effect is of course strongly magnified in electromagnetic processes such 
as prompt photon production because of the valence quark electric charges, as we 
shall see in Section |3l 



2.2.1 Shadowing 

On top of these isospin corrections, parton densities are known to be modified in a 
nuclear environment over the whole Bjorken-x range (see [42-44] for reviews). To 
take into account such shadowing effects, we use the global LO QCD fit of Deep 
Inelastic Scattering (DIS) and Drell-Yan data performed by Eskola, Kolhinen, and 
Salgado (EKS98) [45], who extract the ratio S a /A(x, M) of the nuclear PDF over the 
free proton one: 

C ( AJf\ F a/A(X,M) 

S ^ (X ' M)= AF alr ( X ,MY (2 ' 7) 

The ratio ( |2.7|) depends on each parton species a as well as on the factorization 
scale M through DGLAP evolution. The function S a /A is smaller than 1 at small 
x <C 0.01 (shadowing) and large x > 0.3 (EMC effect), while slightly larger than 
1 (antishadowing) for valence quarks and gluons at x ~ 0.1. Note that the small- 
x region is only poorly constrained by the currently available data, hence quite 
uncertain [46]. 

In the calculations to come, we shall mark the difference between our predictions 
in Au-Au collisions when shadowing effects are taken into account (labelled "Au Au 
EKS98" in the figures) and where they are not ("Au Au"). 

2.2.2 Energy loss 

In the nuclear predictions discussed so far, we assumed implicitly that the hot and 
dense-medium probably produced in high energy heavy-ion collisions does not modify 
either the pion or the prompt photon production process. 

We now suppose that the quarks and gluons produced in the partonic subprocess 
with momentum k ± undergo multiple scattering in the medium. Doing so, they lose 
an amount of energy e with a probability V(e, k ± ). Using a Poisson approximation 
for the soft gluon emission process, this probability distribution is related to the 
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medium- induced gluon spectrum dl/du in [47] and later determined explicitly in 
Refs. [48,49]. In order to make the connection between the energy loss process and 
the quenching finally observed, we follow here the approach of Ref. [27], in which the 
fragmentation variable is shifted from z to z* = z/(l — e/k ± ). The medium- modified 
fragmentation functions in that model thus read [27] 

M r ,k ± ) = r ±(1 de V d (e,kJ z; n D^ /d (z; n} M F ). (2.8) 
Jo 

Let us mention that another attempt at modelling fragmentation functions in the 
medium, accounting for all leading and subleading successive parton branchings, has 
been performed recently within the modified leading-logarithmic approximation [37] . 

The typical amount of energy loss depends on the one scale entering the medium- 
induced gluon spectrum, uj c = 1/2 q L 2 ; here the transport coefficient q is defined as 
the typical kick in transverse momentum space per unit length that the hard gluons 
undergo, and L is the medium path length. While L should be integrated over the 
whole production volume, we shall take here for simplicity a mean length (L) = 5 fm 
in the following calculations. Such a model was shown to describe successfully the 
observed hadron attenuation measured by EMC [50] and HERMES [51] in semi- 
inclusive DIS on nuclear targets [52]. 

Although perturbative estimates for the transport coefficient in a hot quark- 
gluon plasma have been suggested [39], the precise value of u c for the medium pro- 
duced at RHIC energy is somewhat difficult to estimate from first principles. Based 
on the quenching of single-inclusive pion data measured by PHENIX in central Au- 



Au collisions (discussed in Section |3.3|) , we shall use the u c = 20-25 GeV range in 



our calculations. This estimate will be critically discussed in Section [3.4.1 . 
2.3 Ingredients 

The proton parton distribution functions F it j/ p are taken from the NLO (LO) CTEQ6M 
(CTEQ6L) parametrization [53]. The Kniehl-Kramer-Potter (KKP) (N)LO frag- 
mentation functions into neutral pions [54] were used. The Bourhis-Fontannaz- 
Guillet-Werlen (BFGW) [55,56] NLO fragmentation functions for photons were cho- 
sen for both LO and NLO computation in view of the lack of recent leading-order 
determinations of photon fragmentation functions. All scales were taken to be equal, 
H — M — Mp, and chosen so as to minimize the scale-dependence of the NLO 
predictions. In order to investigate the scale sensitivity of the single-inclusive QCD 
calculations, hence part of the theoretical uncertainty 2 , all scales were allowed to 
vary simultaneously in a given range around the optimal scale. The scale-fixing pro- 



cedure in single-inclusive particle production will be discussed in Section [34]. In the 
double-inclusive 7-71" channel, all scales are taken to be given by half the prompt 
photon momentum. 

important uncertainties also arise from the poorly known fragmentation functions. 
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3. Single inclusive pion and photon production 



3.1 Scale-fixing procedure 



Before comparing NLO pQCD calculations with RHIC data, all scales - which should 
be O (p ± ) - have to be fixed within a given prescription. Let p opt be the optimal scale, 
defined as the one which minimizes the scale-dependence of the NLO predictions: 



d da 



dp dp A 



(3.1) 



Mopt 



Although p opt /p ± may in principle depend on p ± , it is fixed here in an arbitrary (but 
sufficiently hard) p ± = 20 GeV value. 




|i / p T n / p T 

Figure 1: LO (dashed) and NLO (solid) predictions for pions (left) and prompt photons 
(right) at p ± = 20 GeV are shown as a function of the renormalization scale, p (normalized 
to p ± ). The factorization scale M and the fragmentation scale Mp are taken to be equal 
to p. The optimal scale is chosen so as to minimize the scale dependence of the NLO 
predictions (see text). 



In Fig. |T] it is shown the scale variation of the LO (dashed) and NLO (solid) pre- 
dictions for pions (left) and prompt photons (right). As expected, NLO predictions 
prove to be much more stable than the LO calculations. In particular, a stable point 
can be found at NLO accuracy while the LO cross sections monotonically decrease 
with the scale p. From Fig. [I], the optimal scale defined in Eq. ([3.1|) proves to be actu- 
ally well below p ± : roughly p? t ~ p ± /2 for pions and p? t = p ± /4 for photons. Since 
we want to perform perturbative calculations down to pretty low p ± ~ 4 GeV values 
in order to compare with data, we shall take slightly larger scales, p* = p ± /\^2 
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and [i 1 — p ± /2, than the strict requirement Eq. ( |3.1|) . It appears moreover that us- 
ing n2 pt = P ± /2 allows for an excellent description of prompt-photon measurements, 
from fixed-target to collider experiments [57]. The theoretical uncertainties discussed 
in the following Section are estimated from the scale-variation around these central 
values. 

3.2 Proton-proton collisions 

The single-inclusive NLO production in p-p collisions is presented in this Section. In 
Fig. ^|are shown the pion (left) and the prompt photon (right) p ± spectra. Predictions 
are given from p ± > 4 GeV - below which perturbative calculations are not expected 
to be reliable - and up to p ± = 30 GeV, above which cross sections are too small to 
be measured with the current RHIC luminosity. At small p ± , pion production cross 
sections prove almost two orders of magnitude, that is O (a s /a), larger than prompt 
photon cross sections. At larger p x > 20 GeV, however, both processes turn out to 
have a similar yield, since the fragmentation mechanism in the pion channel becomes 
kinematically disfavoured with respect to the direct photon process. 




p T (GeV) p T (GeV) 



Figure 2: Single-pion (left) and single-photon (right) invariant cross section at mid- 
rapidity in p-p collisions computed at NLO accuracy, varying simultaneously the factor- 
ization, the renormalization and the fragmentation scales from p a t /\/2 to y/2fi opt . The 
PHENIX data for pions [58] (9.6% normalization error not shown) and the PHENIX pre- 
liminary data for photons [59] are also shown for comparison. Photons and pions are 
produced in the [—0.35; 0.35] rapidity interval. 

The band in Fig. ^| represents the theoretical systematic error of the NLO calcu- 
lation coming from the scale-fixing procedure. In order to discuss this uncertainty on 
a more quantitative level, Fig. ^ displays the ratio of the NLO cross sections at scales 
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Figure 3: Scale dependence of the NLO predictions for single-pion (left) and single-photon 
(right) production at RHIC. Predictions varying simultaneously all scales from /U opt /\/2 to 
v2Mopt are normalized to the "central" prediction n = M = M F = \i . . 
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Figure 4: PHENIX data for single-pion (left) and PHENIX preliminary single-photon data 
(right) normalized to the NLO predictions. The experimental error bars are shown as solid 
(statistical) and dashed (statistical plus systematic) lines, and the theoretical systematic 
uncertainty is shown as a box. The band indicates the 9.6% normalization error in the 
PHENIX measurements. 



varying from /i opt /v^2 to y/2[i opt , normalized to the "central" prediction, all scales 
being equal to // opt (see Sect. |3T| ). This uncertainty turns out to be rather stable as 
a function of p ± for both processes, although the magnitude in both channels differs 
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somewhat. For pions, varying the scales affects cross sections by roughly 30%, while 
a smaller 15% dependence is observed in prompt photon production. The reason is 
twofold: on the one hand, pion production requires an additional a s coupling (hence 
a larger renormalization scale dependence); on the other hand, the direct photon pro- 
cess does not show any dependence on the fragmentation scale. Note that varying 
scales independently (e.g. crossing the scales /i = \/2 p, and M F = A* opt / a/2) leads 
to predictions that lie within our uncertainty band 3 . We should repeat, however, 
that the photon fragmentation functions are much less constrained than those in the 
pion sector; this leads to an additional theoretical uncertainty, not quantified here. 

Let us now discuss the pion and photon NLO predictions in comparison with the 
recent PHENIX measurements [58-60] . As can already be seen in Fig. [|, an excellent 
agreement is found on the whole p ± range. To be more precise, the data over the 
NLO theory ratio is computed in Fig. |j. Boxes indicate the theoretical uncertainty 
on this ratio coming from the scale dependence, and the vertical error bars show the 
statistical uncertainty in PHENIX data, with (dashed) or without (solid) systematic 
error. 

For pions, it is worth noting that errors coming from theory and experiment 
remain not too large over the whole p ± ^ range. At small p ± ^ , the main uncertainty is 
given by the scale dependence of the NLO computation (approximately 30%) while 
experimental error bars prove as small as 10%. On the contrary, at large p ± ^ , the 
error on the data-over-theory ratio is mostly statistical, because of the low counting 
rates with the present RHIC luminosity. Fig. |] clearly indicates that the present 
NLO 7T° cross sections overestimate the PHENIX measurements by roughly 30%; 
yet the data-over-theory ratio proves to be flat over the whole p ± range. Slightly 
larger scales, say all equal to p Xjr , would significantly improve the description of 
PHENIX data. It is nevertheless quite remarkable to observe such a fair agreement 
between the NLO predictions and the data since the KKP fragmentation functions 
used here have not been constrained by any hadron-hadron scattering data [54]. It 
may therefore certainly be useful to reverse the logics and to use the neutral-pion 
PHENIX data in a global (N)LO QCD fit analysis to further constrain the pion 
fragmentation functions and thus gain additional accuracy on the LHC predictions. 

Let us move to prompt photon production (Fig. |], right). The data-over-theory 
ratio is now completely dominated by the statistical uncertainty in the PHENIX 
data at all transverse momenta. Theoretical predictions prove to be in excellent 
agreement with data 4 . We hope that the experimental errors can be reduced in the 
future, so as to test more drastically the current NLO prompt photon predictions, 
whose fragmentation functions are not well constrained yet. 

3 P. Aurenche and M. Werlen, private communication. 

4 On the contrary, the fj, = M = M F = p ± prescription shown in [59] somehow underestimates 
the PHENIX measurements. 
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3.3 Nucleus— nucleus collisions 

We now come to the production of single- 7r° and single-7 production in Au-Au col- 
lisions. In Fig. |5] is plotted the expected quenching of these two probes (respectively 
on the left and on the right panel) assuming (i) no nuclear effect to be at work 
(labelled "Au Au", solid), (ii) including shadowing corrections ("Au Au EKS", dash- 
dotted), and (iii) including both shadowing and energy loss processes ("Au Au EKS 
20 < uj c < 25 GeV", band). 

While isospin effects in the pion channel are negligible at RHIC, we note that 
antishadowing tends to enhance the nuclear production by roughly 20% in the lowest 
p ±7r bin. Conversely, at large p ±7v > 15 GeV, the large-x EMC effect slightly quenches 
the pion yield. The single-pion spectrum in Au-Au collisions including shadowing is 
therefore somewhat softened with respect to the p-p scattering case, although the net 
effect remains small. In contrast, the effect of the parton energy loss process, shown 
as a band, is dramatic. The quenching factor starts around 0.1 at p^ = 4 GeV 
and smoothly increases to 0.3 at p^ = 20 GeV, in rather good agreement with 
the PHENIX preliminary data [61]. Note that, even if the p ± ^ dependence of the 7r° 
quenching is somehow flattened by the shadowing contribution, no hint for such an 
increase is seen in the data. The rather flat behaviour observed there may be due to 
the geometrical bias introduced when integrating over all possible path lengths [34], 
not performed here. 

In contradistinction to the pion case, the isospin correction is significant in 
prompt photon production (right, solid). Photon production being an electromag- 
netic process, cross sections depend on the light-quark electric charges and are thus 
disfavoured in a nucleus target, less rich in up quarks than a proton is. It is worth 
stressing that the quenching - without any shadowing nor energy loss effects - proves 
as large as 20% at p ± ^ = 20 GeV! Antishadowing slightly compensates the isospin ef- 
fect to make the quenching factor closer to unity below 10 GeV. Unlike the pion 
channel, antishadowing extends up to p ± ^ ~ 20 GeV, at which it is negligible. 
This is so because the Bjorken x probed in the direct photon channel is smaller, 
x — O \ 2p x ^/y/sji than in the pion fragmentation process, x — O (2p Xjr / 'z^yfs). 

Let us now comment on the energy loss effects. Because of the dominance of the 
direct process unaffected by the medium (keeping in mind that its relative contribu- 
tion is scale-dependent and thus somehow arbitrary), the photon quenching is not as 
pronounced as that of the pion. At 4 GeV, prompt photon production is suppressed 
by 40% but remains of order 30% in the largest considered p ± bins. While energy 
loss effects in the present model lead to an increase of the quenching factor with p ± ^ , 
the rather flat behavior seen in Fig. |5] actually comes from the interplay of the energy 
loss process on the one hand and of the isospin effects on the other hand. 

The quenching of single-photon production is unfortunately not yet available due 
to the too low statistics in p-p collisions. In Fig. ^ (right) is shown the ratio of the 
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5 10 15 20 5 10 15 20 



p T (GeV) p T (GeV) 

Figure 5: Ratio of Au-Au over p-p single-pion (left) and single-photon (right) production 
cross section. Calculations are done at LO, assuming (i) isospin (solid), (ii) isospin and 
shadowing (dotted), (iii) isospin, shadowing and energy loss (band) effects (see text for 
details). The PHENIX preliminary data on the single-pion [61] and single-photon [62] 
production in Au-Au collisions normalized respectively to p-p scattering 7r° data and to 
the present NLO prompt photon prediction in p-p collisions are also shown for comparison. 



Au-Au PHENIX photon measurements [62] to the present NLO calculations in p-p 
collisions. In that respect, this is not - strictly speaking - the same quenching factor 
as the one determined in the pion sector for instance. Nevertheless, given the good 
agreement between PHENIX data and NLO calculations discussed in Section |3.2| (cf. 
Fig. ||, right), we believe this ratio to be indicative of the genuine prompt photon 
quenching (i.e. normalized to p-p data). Interestingly, we notice that the photon 
quenching factor turns out to be in very good agreement with the PHENIX measure- 
ments. Although the presently too large error bars do not allow one to disentangle 
the predictions with or without energy loss effects, it is worth stressing that our 
estimated 30% suppression is not inconsistent with these preliminary measurements. 
It would indeed be particularly intriguing not to observe a suppression - though less 
spectacular than for the pions - due to the parton energy loss. Calculations done at 
RHIC using isolation criteria [59] indeed indicate that roughly 20% of the inclusive 
photon yield actually comes with an important hadronic activity 5 ("jet-like pho- 
tons"). Let us mention that our predictions follow the same trend as the calculation 
by Jeon, Jalilian-Marian, and Sarcevic [63] who first attempted to determine prompt 
photon quenching at RHIC. However, we can regret the lack of clear relationship be- 
tween the energy loss probability distribution used in [63] and the medium-induced 

5 I thank P. Aurcnche and S. Bathe for discussions on this issue. 
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K ~f/-K° = Z7ZZ — ~~ Znv — ~Z^T\ • t 6 - 2 ) 



gluon radiation computed in QCD by BDMPS [39]. 

The difference between the single-pion and the single-photon quenching discussed 
above can be better seen in the ratio of the photon total yield - i.e. including prompt 
photons as well as the background photons coming from the decay of neutral pions - 
over that very background: 

a(pp — > 7 X) + a(pp — > 7T° X — > 7 X) 
a(pp — > 7r° X — > 7 X) 

The 7T° decay contribution 7r° — > 7X to the photon total yield may be simply esti- 
mated from the slope n of the single-pion spectra assuming a power law behaviour 6 

da I 'dp ± ~p7 n , [4,64] 

a(pp^n°X^ 1 X) _ 2 
cr (pp — >• 7t° X) n — 1 

In p-p collisions (Fig. || solid line), the ratio increases from 1 at small transverse 
momentum - that is with a negligible prompt photon production - up to 2 at 20 GeV, 
above which the prompt photon signal takes over the background. Shown respectively 
as a dash-dotted and as a dotted line are the same ratios in Au-Au collisions when 
adding isospin effect, respectively with and without shadowing corrections. The 
overall effect of these cold nuclear matter corrections is pretty weak, making the 
ratio slightly smaller than what is expected in p-p collisions. 

On the contrary, one might expect dramatic effects of parton energy loss on 
this observable. Indeed, since the pion yield proves more suppressed by this process 
than does the photon yield, this ratio should increase much faster with p , in the 
presence of the hot medium. The ratio has been computed assuming the previous 
uj c = 20-25 GeV range in Fig. |6| (band). As expected, Rj/^o turns out to be almost 
twice as large at 20 GeV in Au-Au collisions with energy loss than in p-p collisions. 
Moreover, it is particularly interesting to notice that this prediction proves to be in 
good agreement with the recent PHENIX data in central Au-Au collisions, although 
the predicted slope tends to underestimate that of the data. Nevertheless, we find 
it quite remarkable that both single inclusive 7r° and 7 channels can be described 
within the same model and with one common parameter. This allows us to use such a 
model further and to consider, in Section f|, photon-tagged correlations as a possible 
tool to probe more quantitatively the medium produced at RHIC. 

3.4 Discussion 

3.4.1 RHIC energy density 

The good agreement between the single-inclusive p-p and Au-Au scattering data and 
the theoretical calculations is a first basic check of the present approach. Perhaps 

6 The slope n, and therefore the ratio -R decay , is estimated from the PHENIX preliminary data in 
p-p collisions [58]. Since the quenching factor is remarkably flat above 4 GeV, n should not change 
much from p-p to Au-Au collisions. 
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Figure 6: Ratio of the overall photon production over the background photon production 
coming from the 7r° decay. LO calculations are done in p—p and Au-Au collisions (see text 
for details). The PHENIX data [62] are shown for comparison. 



more importantly, this allows us to constrain the unknown parameter of the model, 
u c , and then to relate it to other physical quantities regarding the dense-medium 
produced. We would like here to briefly comment on its absolute value in order 
to get a rough estimate of the energy density currently achieved in central Au-Au 
collisions at RHIC. 

The dynamical scaling law sets the relation between the time-averaged (q) and 
the initial-time q(tn) transport coefficient [49]: 



<$> = J? [ dt(t-t ) (jj q(to) 



) q(t ) when t < L, (3.4) 



2-a \L 

accounting for the expansion of the produced medium with density n(t) oc t~ a . 
Assuming in the following a purely longitudinal expansion (a = 1) in (|3.4j) , the 
transport coefficient at an initial time to is thus given by [65] 

m * ^j- (3.5) 

In Li 



Using the estimate uo c = 20-25 GeV in ( |3.5|) , one then gets the transport coefficient 
q(to = 0.5 fm) ~ 1.6-2 GeV 2 /fm at an early time to — 0.5 fm after the reaction. 

In Ref . [65] , Baier related the transport coefficient computed perturbatively [39] 
to the medium energy density, e, both for hot pion gas or quark-gluon plasma. 
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Fig. [7| exhibits the generic power law dependence of the transport coefficient on e. 
In order to make a crude estimate on the energy density reached in central Au-Au 
collisions at RHIC energy, we superimposed on this curve our estimates for the time- 
averaged (open circle) and initial time (full circle) transport coefficient. This leads 
to an energy density of about e ~ 15 GeV/fm 3 , i.e. 100 times that of cold nuclear 
matter. For comparison, also shown on this plot is the nuclear matter transport 
coefficient, both the perturbative estimate (open square) [39] as well as the slightly 
larger but consistent value extracted from Drell-Yan production in pion-induced 
nuclear collisions (full square) [66]. 
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Figure 7: The solid (respectively, dashed) line represents the energy density dependence 
of the transport coefficient in a thermally equilibrated quark-gluon plasma (respectively, 
pion gas) [65]. The squares indicate the transport coefficient for cold nuclear matter com- 
puted perturbatively (open square) [39] and extracted from Drell-Yan and DIS data (full 
square) [52,66]. The open and full circles represent respectively the mean (q) and the 
initial-time q(to = 0.5 fm) transport coefficient at RHIC (figure adapted from [65], with 
permission). 

Let us recall here that our estimate is no more than a guess, as it may somehow 
depend on our specific assumptions, regarding for instance the longitudinal expan- 
sion, the choice to take a fixed medium length, not to say the intrinsic uncertainty 
coming from the perturbative calculation of the transport coefficient. Nevertheless, 
we notice that our result, e HHIC /e cold — 100, proves to be on the same order of magni- 
tude as estimates based on soft particle production at mid-rapidity [67]. Moreover, 
this number should be seen as a lower estimate, since our assumptions - in partic- 



- 15 - 



ular the geometrical modelling as well as the absence of a vanishing cut-off in the 
BDMPS medium-induced gluon spectrum - is likely to underestimate the initial-time 
transport coefficient q(t ). 

It would also be interesting to get an estimate on both temperature and energy- 
density from the above transport coefficient within a quasi-particle picture. The 
gluon damping rate 7, or equivalently its inverse mean free path A 9 _1 , has been 
computed in SU(iV c ) colour gauge theory in the weak coupling limit [68]. It is given 



Fitting the entropy computed on the lattice in quenched QCD, Peshier extracted the 
coefficient appearing in the logarithm of (|3.6|) , c = 2.23 [69]. From the leading-order 
Debye mass m D = y/4ira s T, we obtain the relation between transport coefficient 
and temperature: 



Taking the transport coefficient q = 1.6-2 GeV 2 /fm and a s = 1/2 at rather soft 
scales, the temperature of the medium is T = 355-385 MeV. In a gas of weakly 
interacting gluons, this would correspond to an energy density e RHIC = 16 n 2 T 4 /30 = 
11-15 GeV/fm 3 , which compares well with the former estimate. 

3.4.2 Uncertainties on medium-induced photon production 

Despite the good agreement between the model predictions and the PHENIX data, 
the currently large experimental and theoretical uncertainties do not rule out other 
possible channels for photon production in a hot medium. 

To leading order, the photon can participate directly in the hard partonic subpro- 
cess and thus escape without any strong interaction. On top of that direct channel, 
photon production may also come from the collinear fragmentation of quarks and 
gluons produced in the hard process. Although of higher order according to power 
counting rules, the large logarithm ln(Q 2 /A 2 ) ~ a' 1 coming from the integral over 
the photon transverse momentum makes this process contribute to the same leading 
order O (a a s ) in the perturbative expansion as the direct process [11]. The formation 
time needed to produce such a photon-parton system with a small invariant mass 
exceeds by far the typical lifetime of the hot medium. This allowed us to consider 
a two-stage process: the multiple scattering of the hard parton in the medium is 
followed on a much larger time-scale by the parton-to-photon fragmentation process 
in the vacuum. 

However, the possibility for the photon to be produced while still inside the 
medium has not been considered here. In particular, the multiple scattering in- 
curred by hard partons induces the emission of soft gluons as well as soft photons. 



by 




(3.6) 
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Within his path-integral picture for the parton energy loss mechanism [70,71], Za- 
kharov recently computed this medium-induced photon bremsstrahlung contribution 
at RHIC [72]. The enhancement of photon production is particularly noticeable in 
the moderate p ± _^ range, say below 20 GeV, and could therefore somewhat balance 
the photon quenching predicted in Fig. [5|. In addition, large transverse momentum 
partons may couple to the thermal quarks and gluons in the medium through Comp- 
ton scattering or qq annihilation. This jet-photon conversion mechanism, considered 
by Fries, Miiller and Srivastava in Ref. [73], is particularly rich since it would di- 
recty reflect the hard parton instead of the measured hadron spectrum, provided this 
channel is dominant in photon production. 

Several analyses came out, recently, which embed some - or all - of these channels 
into hydrodynamical evolution [36,74-76]. The agreement with the recent PHENIX 
data [59] is remarkable down to low p ± values, p ± > 1 GeV. In particular, it appears 
that the new preliminary and more precise data [6] are better reproduced with a 
thermal contribution [74]. This is clearly an important observation, which would 
deserve further study. Note, in passing, that a strong depletion of the perturbative 
QCD yield expected at low p ± ^ because of the energy loss mechanism would make 
the photon enhancement even much stronger! However, one should not expect the 
NLO pQCD calculations to be reliable and well constrained at such a small p ± . In 
order to bypass the comparison with the pQCD calculation extrapolation, it would 
be extremely interesting to measure the quenching factor down to 1 GeV and to see 
whether the photon enhancement remains. 

Each of these individual mechanisms is unfortunately poorly under control, not to 
mention the uncertainty due to the space-time evolution of the expanding medium. 
Given this variety of processes, it appears that the medium effects on (prompt) 
photon production still remain quite unknown and model-dependent. Hopefully more 
precise data in the near future will shed light on the prompt photon production 
mechanism in nucleus-nucleus collisions. In the following, however, we shall discuss 
double inclusive pion-photon correlations at some specific regions of phase space in 
which the photon is (mostly) produced directly in the hard process. Consequently, 
our predictions should fortunately not depend much on the medium-induced photon 
production process. 



4. Double inclusive photon-pion production 

The momentum imbalance 



spectrum between a hard prompt photon and a much softer (but still hard) hadron 
produced in hadronic collisions may allow for the determination of the hadron frag- 
mentation function, D^(z ~ z 34 ), the photon transverse momentum balancing that 
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of the parton i, which fragments into the hadron. At least this two-body kinematics 
may be a valid picture when higher-order corrections, briefly discussed below, re- 
main small. Moreover, the fixed-order calculations should not be reliable at large z 
because of the large logarithms a s n m 2n (l — z j7r ) and a s n ln 2n_1 (l — z 17T ), due to soft 
and collinear gluon emissions, which need to be resummed to all orders. 

This triggered a recent phenomenological study of various pion-photon and 
photon-photon correlations in heavy- ion collisions at LHC energy [28] . Our goal here 
is to investigate similarly such momentum correlations at RHIC energy, restricting 
ourselves to the 7-71" channel. For more details concerning the other kinematical 
variables used here, see [28]. 

4.1 Proton-proton collisions 

In order to probe the pion fragmentation function efficiently through the 7-71° 
momentum imbalance spectrum, asymmetric cuts p cut <C p cnt are required to make 
the range covered in the z fragmentation variable as wide as possible. In addition, 
the pion momentum needs to be hard enough to ensure the perturbative regime to be 
at work, while the photon momentum should not be too large to maintain reasonable 
counting rates: 

a qcd « P ^ « « v'i /2. (4.i) 

In the following we choose the cuts p cnt = 3 GeV and p cnt = 10 GeV. 

J-7T -L-7 

Fig. H shows the pion and the photon p ± spectra (top) as well as the invariant 
mass m 77r and the transverse momentum of the pair q ± ^ (bottom) computed in p-p 
collisions at LO (open circles) and NLO (full squares) accuracy. To leading order, the 
two particles are emitted back to back in the parton-parton center-of-mass frame; 
hence the invariant mass distribution has a threshold given by 2 /pcutpcut ~ n GeV 

and the q ± ^ spectrum shows a maximum around p™* — p° ut = 7 GeV. At NLO, 
new configurations in momentum space show up since the pion and the photon may 
be emitted at a smaller relative azimuthal angle (we apply a minimal cut 77r > 
7r/2). This fills in particular the small invariant mass region (with a new threshold 
/pcutpcut r-j 5 5 GeV) and shifts the q T spectrum to slightly larger values. Except 

y J-7T -L'y 1 

in this region of phase space, however, LO predictions appear quite stable with 
respect to the small NLO corrections, say roughly 20-30% at most. We should 
however keep in mind that the strength of NLO corrections strongly depends on 
the renormalization, factorization and fragmentation scales assumed in the pQCD 
calculation (here, all fixed at p ± /2). 

The distribution in the 7-71"° momentum imbalance z is also determined to 

' 7"7T 

LO and NLO accuracy (Fig. |]). As already pointed out in Ref. [28], its shape turns 
out to be reminiscent of the pion fragmentation function. Above the ratio of the 
cuts, p cnt /p cnt = 0.3, the distribution decreases as the pion momentum and thus the 
fragmentation variable z gets larger. Similarly, the distribution is strongly suppressed 
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Figure 8: The pion p ± ^ and the photon p ± ^ transverse momentum distributions (top), 
the invariant mass m^, and the pair transverse momentum q ± ^ (bottom) are com- 
puted in proton-proton scattering to LO (open circles) and NLO (full squares) at = 

200 GeV. Both the photon and the pion are produced at rapidity [—0.5; 0.5] and the cuts 
pcut _ -j^q Q e y anc [ pcut _ 3 Q e y are imposed. 



below z in < 0.3, since larger photon momenta are needed, keeping p ± ^ close to its 
lower cut. Again, higher-order corrections prove large when the two particles are 
emitted at 7ir > ir/2, thus at small values of z^ (recall that z i7i oc cos0 77r ), but 
moderate elsewhere. Let us furthermore insist that, despite the regular behaviour 
of the NLO predictions, soft gluon resummation may affect the distribution at high 



4.2 Nucleus— nucleus collisions 



In Fig. [10| the quenching factors of the p^ , p ± ^ , q ± ^ and LO spectra in Au-Au 
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Figure 9: The 7-71"° momentum imbalance z 7Jr distribution is computed in proton-proton 
scattering to LO (open circles) and NLO (full squares) at yfs = 200 GeV. Both the photon 
and the pion are produced at rapidity [—0.5, 0.51 and the cuts » cut = 10 GeV and p cut = 
3 GeV are imposed. 

collisions are computed, including nuclear shadowing and with energy loss (u c = 
20 GeV, open squares) or without it (u c = GeV, full squares). The effect of isospin 
and nuclear shadowing is rather small for all observables when both the pion p ± ^ and 
the photon p^ momenta are close to their respective cuts, i.e. where distributions 
are maximal. However, the quenching factor decreases down to 0.7 when momenta 
become larger (e.g. at p ^ ~ 40 GeV) since higher x are probed in the Au nuclei. Just 
as the single production case, the quenching comes from the interplay between the 
lack of up quarks in the nuclei and the nuclear EMC effect. Note that the depletion 
is also pronounced at small q ± , when the pion momentum is of the order of the 
photon momentum, p ± ^ ~ p ± ^ . In this specific region, not only the pion but also the 
photon are produced by collinear fragmentation. This double fragmentation process 
then requires highly energetic (i.e. large-x) partons to be produced to fragment into 
the two detected particles, therefore leading to a similar suppression. 

The effect of the energy loss on each of the spectra is more dramatic. In particu- 
lar, unlike isospin and shadowing, the smaller the momenta the stronger the effects. 
As an example, the quenching is as large as 0.3 close to the photon momentum cut, 
p ± > 10 GeV, or even 0.2 at small pair momentum, q x ^ ~ 5 GeV. As in the single 
production case, energy loss effects naturally die out at larger momenta, since the rel- 
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Figure 10: Same distributions as Fig. 8. Distributions are computed in Au-Au collisions 
at y/s = 200 GeV assuming shadowing (full squares) and shadowing with energy loss (open 
squares) normalized to the p-p distributions. 



ative energy loss, e/k ± , gets smaller. Once more, though, the competition between 
energy loss on the one hand and isospin/shadowing on the other leads to various 
patterns for the quenching. While the invariant mass quenching remains remarkably 
flat on the whole kinematic range, the quenching slightly decreases (respectively, in- 
creases) with p ±jr (respectively, p ± ). Particularly interesting is the q ± ^ spectrum 
behaviour in the medium. As stressed above, increasing above its maximum (at 
p cut — p cut ) amounts to increasing the photon momentum. This is the reason why 

J-7T J-7 

both the p ± ^ and the spectra are similarly quenched (Fig. 0, left). At very 
small and negative q ±i7v , however, the pion becomes harder while the photon mo- 
mentum remains of order p cut . Since the double fragmentation process makes energy 
loss effects weaker, the quenching turns out quite naturally to be maximal around 
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Let us finally discuss the 7-71 momentum imbalance quenching shown in Fig. 1 1 



Remarkably, the effect of isosp in/shadowing is completely negligible (much less than 
10%) above the cut ratio, z > p cut /p cnt . In that domain, the pion momentum 
p ±w grows while the photon momentum - and therefore the transverse momentum 
of the parton that fragments, to leading order - is kept fixed. The typical values 
of x do not change much and hence neither does the quenching due to isospin or 
shadowing. These effects are then located in the small region z < p cnt lp cnt where 
the photon momentum is much larger, p ± ^ ^> p^ ut . The energy loss mechanism 
leads to a completely different picture for the quenching. At larger z ~ z, the 
phase space for gluon emission is dramatically restricted (see Eq. Q2.8|) ) and the 
quenching becomes more pronounced. This explains the decreasing behaviour of the 
quenching factor, from 0.5 at z j7T ~ 0.1 down to 0.15 at z ~ 0.8. At very large 
z iiT > 0.8, the quenching factor starts increasing again because of the onset of the 
double fragmentation process, as already mentioned. 

It is clearly the momentum imbalance variable that offers the largest observable 
difference whether including energy loss (u c = 20 GeV) or not (u c = GeV) in the 



model. Fig. |TT] also justifies a posteriori our choice of extremely asymmetric cuts 
between the pion and the photon momenta, in order to isolate as much as possible 
the effects of isospin and nuclear shadowing. 
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Figure 11: The 7-71" momentum imbalance z 11T distributions is computed in Au-Au 
collisions at ^/s = 200 GeV assuming shadowing (full squares) and shadowing with energy 
loss (open squares) normalized to the p-p distributions. 
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4.3 Counting rates 

The momentum correlations between a (ideally directly produced) prompt photon 
and a pion, as discussed above, are clearly of interest at RHIC energy in order to 
probe and constrain the present energy loss models, which have described successfully 
the single pion/photon production. Hopefully double inclusive 7-71- data can be 
available in the near future at RHIC. 

For this reason, we would like to discuss here what the expected counting rates 
are for such processes. The number of events is given by [4] 



a 



gco 



\ /"hard r I at \\ AuAu _hard n 
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The RHIC highest integrated luminosity delivered in Run 4 and measured by PHENIX 
is taken to be C = 0.7 nb _1 [77]. Cross sections in Au-Au collisions range roughly 
from 10~ 4 to lCT 1 nb/GeV, so that approximately M = 2-2 10 3 events/GeV are 
expected. Regarding the imbalance distribution, rates should be dJ\f/dz = 5 10 2 - 
5 10 4 events. These numbers are given obviously without considering any additional 
kinematic cuts or acceptance restrictions. A higher luminosity at RHIC is clearly 
hoped for, as it would certainly allow for a more systematic investigation of 7-71- or 
7 — 7 correlations with a great variety of cuts. 



5. Summary 

Single inclusive production of pions and photons at large transverse momentum 
ip± > 4 GeV) has been computed perturbatively in p-p and Au-Au collisions at 
RHIC energy (y^i = 200 GeV). In p-p collisions, NLO calculations compare success- 
fully with the PHENIX data on the whole p ± range for both pion and photon spectra. 
Theoretical uncertainty from the scale fixing in the calculation was also examined. In 
Au-Au collisions, we determined the p ± spectra, assuming possible shadowing cor- 
rections in the nuclear parton densities or energy loss effects on the fragmentation 
process. The quenching of single-pion production observed by the PHENIX collabo- 
ration can be reproduced, assuming the typical scale u c = 20-25 GeV for the energy 
loss process. From these values, we estimate the energy density reached at RHIC 
in central Au-Au collisions to be e RHIC > 10 GeV/fm 3 at an initial (and somewhat 
arbitrary) time to — 0.5 fm. Within the same energy loss model, the quenching factor 
proves much less pronounced in the photon sector, because of the presence of the 
direct process channel unaffected by the medium. The expected photon quenching is 
found to be roughly 20%, therefore in agreement with the PHENIX measurements. 
Finally, the ratio of the photon total yield over the pion decay background was found 
to fairly reproduce the present data. 
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The production of double inclusive 7-71" momentum correlations has then been 
investigated in detail. In p-p scattering, LO predictions were shown to be quite 
stable with respect to higher-order corrections, except in some specific regions of 
phase space. It is the case in particular when pions and photons are emitted at small 
relative azimuthal angle, thus at small momentum imbalance. Using the asymmetric 
cuts p cut = 3 GeV and p cut = 10 GeV, the quenching of various correlation spectra in 

J-7T -L7 

Au-Au collisions are predicted. The momentum imbalance distribution z iiT is seen 
to be particularly sensitive to the medium-modified fragmentation dynamics. The 
expected counting rates assuming the largest RHIC integrated luminosity in Run-4 
are encouraging, even though a higher luminosity would be required for a thorough 
study of photon-tagged correlations at RHIC. 

As an outlook, it would be interesting to go beyond the present calculation and to 
perform a systematic comparison of photon-tagged versus hadron-tagged momentum 
and azimuthal correlations. Since the latter is shown to be quite sensitive to the 
surface emission, and hence to the geometry of the collision [34], one could naturally 
expect significant differences in such correlations when triggering on prompt photons. 
Going from a coloured (parton or fragmentation photon) to a blind (direct photon) 
trigger would scan different energy density profiles. In that respect, performing 
isolated versus non-isolated photon correlation measurements would be ideal. 
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